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INTRODUCTION 

There  has  been  an  increasing  interest  in  the  use  of  mode  tuned  or  stirred 
reverberation  chambers  for  conducting  electromagnetic  coupling  and  upset 
experiments.  Mode-stir  chambers  consist  of  a  high  quality  factor  (Q) 
metallic  reverberation  chamber  in  which  a  stirring  device  such  as  a  metallic 
paddle  wheel  is  used  to  incrementally  or  continuously  alter  the  boundary 
conditions  within  the  chamber.  The  goal  is  to  achieve  an  isotropic 
homogeneous  field  everywhere  in  the  volvuse  of  the  cavity  except  near  the 
walls.  Field  hlomogeneity  is  achieved  in  this  manner  by  exploiting  the 
pseudo-statistical  nature  of  each  of  the  eigenfxinctlon' s  contribution  to  the 
field  level  at  a  given  point  well  within  the  chamber  volume.  Extensive  work 
has  been  done  by  the  National  Institute  of  Standards  and  Technology  (NIST) 
(previously  the  National  Bureau  of  Standards)  and  the  Naval  Surface  Weapons 
Center  to  both  optimize  chamber  design  and  characterize  chamber  fields.  The 
NIST  has  achieved  field  uniformities  of  <  ±  3  dB  above  1  GHz  and  <  ±  2  dB 
above  2  GHz  in  a  2.74  x  3.05  x  4.73  m  welded  steel  enclosure  using  the  paddle 
wheel  technique  [Ref.  1] . 

Reverberation  chamber  testing  offers  several  advantages  over  aneehoic  or 

plane  wave  illumination  testing.  For  example,  E-fleld  levels  in  the 

.  /  • 

thousands  of  volts  per  meter  can  be  achieved  using  modestly  low  power  sources 
such  as  200  U  TWT  amplifiers.  Coupling  cross  sections  can  be  measured 
independent  of  the  angle  of  incidence  (this  is  particularly  beneficial  when 
testing  subsystems  that  would  normally  exi$t  inside  an  equipment  bay  or 
cavity) .  On  the  other  hand,  dependence  on  mechanical  mode  stirring  can 
complicate  data  acquisition  and  interpretation  as  well  as  require  long 
periods  of  time  (as  much  as  10  hours  per  test)  to  acquire  a  complete  data 
set.  These  complications  are  primarily  due  to  the  requirement  that  data  be 
measured  for  many  different  positions  of  the  paddle  wheel  for  each  frequency 
and  power  level  because  field  uniformity  inside  the  volume  of  the  chamber  can 
only  be  obtained  by  averaging  the  fields  Influenced  by  each  position  of  the 
tuner  over  many  positions.  Hence,  real  time  field  uniformity  can  only  exist 
on  time  scales  oj.the  same  order  as  the  paddle  wheel  rotation  rate. 
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Flgtire  1.  The  NT  Enclosure. 


Figure  2.  The  GIFT  box  enclosure. 


Ic  was  the  purpose  of  this  experiment  to  investigate  the  use  of  electronic 
mode  stirring  in  a  reverberation  chamber  to  perform  both  coupling  and  upset 
testing.  While  the  conventional  method  of  mode  stirring  holds  the  frequency 
constant  and  varies  the  cavity's  boundary  conditions  In  order  to  obtain  a 
sufficiently  large  sample  of  eigenfunction  contributions  to  the  field  levels 
at  a  given  point  In  the  chamber,  the  method  described. here  does  not  vary  the 
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boundary  conditions,  but  instead  averages  the  eigenfunction's  contributions 
o’fer  a  narrow  band  of  frequency.  Expected  advantages  of  this  technique  over 
the  more  conventional  method  included  shorter  test  times  (from  hours  to  less 
than  5  minutes),  simplified  data  acquisition  and  control,  and  more 
interpretable  results.  Because  of  their  availability  and  in  order  to 
determine  the  effect  chamber  size  has  on  Che  technique,  two  separate 
reverberation  chambers  were  chosen,  one  relatively  small  and  one  relatively 
large  (Figures  1  and  2) .  The  larger  enclosure  is  referred  to  as  the  NT  ■ 
enclosure  because  it  was  borrowed  from  the  Nuclear  Technology  (NT) 
Directorate  here  at  the  Weapons  Laboratory  (now  the  Phillip.s  Laboratory) , 
while  the  other  is  referred  to  as  the  GIFT  box  (an  acronym  from  a  previous 
test). 
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THEORETICAL  DEVELOPMENT 

B-DOT  PROBES 

There  are  several  different  sensors  available  for  measuring  electric  and 
magnetic  fields.  B*dot  probes  were  chosen  for  this  experiment  because  of 
their  availability,  small  size,  and  broad  bandwidth.  The  voltage  output  of 
the  B-dot  probe  is  directly  proportional  to  both  the  probe's  effective  area 
A,,  and  the  time  rate  of  change  of  the  B  field  dotted  into  the  normal  vector 
to  the  plane  of  the  probe's  loop; 


v(t) 


(1) 


Since  only  continvious  wave  signals  are  used  in  this  experiment,  it  is  more 
useful  to  consider  the  Fourier  transform  of  the  above  expression: 

G(i«)  •  r  ■ 

rf  •  B(Jti) 

where  G(Ju)  is  the  transfer  function  that  relates  the  voltage  output  to  the 
amplitude  and  phase  of  the  B  field.  At  a  fix.«d  frequency,  w.,  the  power 
density,  P^,  due  to  the  component  of  the  i  field  normal  to  the  probe  loop,  B, 
can  be  related  to  the  probe's  power  output,  P.  into  a  50  0  instrument  as 
follows : 


Since  ■ 


1  egg 

2  Ho 


2 


1 

SO 


50c 


(3) 


where  c  la  the  speed  of  light  and  ft,,  is  the  permeability  of  free  space  (the 
medium  is  assumed  to  be  a  vacuum) .  The  expression  for  the  power  density,  P„, 
above  was  found  by  multiplying  the  energy  density  of  a  plane  wave  by  the 
speed  of  light. 

RECTANGUI/il  CAVITIES 

Electromagnetic  theory  as  applied  to  conductive  rectangular  cavities  can  be 
found  in  many  standard  electromagnetic  texts  and  reports.  For  convenience. 
Several  of  the  more  important  results  are  repeated  here. 
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The  frequencies  at  which  eigeranodes  can  'xist  in  a  rectangular  cavity  of 
dimensions  a,  b,  and  d  are  [Ref.  2,  pp  345]: 


<4) 


where 

m  n,l  -  Any  set  of  integers  >0  except  {0,0,0). 
c  —  Speed  of  light. 

-  Relative  permeability  and  permittivity  of  the 
medium  filling  the  cavity. 

Of  more  interest,  however,  is  the  total  number  of  eigeracodes  that  can  exist 
in  a  small  interval  of  frequency  with  a  center  frequency  of  f.  This  quantity 
is  referred  to  as  the  specific  mode  density.  Blackbody  mode  density 
calculation  results  are  directly  analogous  to  this  problem  and  express  the 
mode  density  a.s  [Ref.  3): 


cW  -  df  (5) 

c’ 

Mode  density  is  important  because  it  determines  the  effective  sample  size  of 
the  pseudo-random  eigenfunction  amplitudes  at  the  point  of  interest.  The 
larger  the  effective  sample  size,  the  better  the  statistical  randomness  of 
the  eigenfunction  amplitudes  at  a  given  point  and  the  better  the  field 
homogeneity.  The  theoretical  mode  density  for  a  1-m’  cavity  with  several 
different  frequency  agile  source  bandwidths  is  shown  ir>  Figure  3. 

Another  important  parameter  is  the  quality  factor  Q: 

o  .  ■«-  Energy  Stored  ^  Energy  Stored  /g\ 

*  Power  Lose  /  Hertz  ^  Pj. 

where  Pj  is  the  net  power  being  transmitted  into  the  enclosure  and  Wg  is  the 
radian  frequency  of  the  exciting  source. 

'S'.',  ■ 
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MODES  IN  A  UNIT  CAVITY 


Centex  Pze^uency  (GHz) 


Figure  3.  Mode  density  as  a  function  of  the  agility  bandwidth 
and  center  frequency  for  a  l-n>®  cavity. 


Composite  0  of  a  Rectangular  Cavity 


The  expression  given  above  fot  the  Q  of  a  cavity  is  singular  in  frequency. 

In  other  words,  the  equation  assumes  that  is  a  monochromatic  source  with 

frequency  Wq.  If  instead  the  spectrum  of  Pj  is  spread  across  a  narrow  band 
of  frequency  (as  is  the  case  for  this  experiment) ,  it  is  convenient  to  derive 
an  expression  for  a  composite  Q  that  jpresents  the  net  effect  of  having  many 
eigenmodes  stimulated  simultaneously.  More  specifically: 


ay 
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where 


Q»<iv 

Wi.Pii 

?TT 

Wt 

"e 


Composite  Q. 

Energy  in  and  power  lost  to  eigenmode  i. 
Total  net  power  transmitted  into  the. cavity - 
Total  energy  stored  in  all  eigenmodes. 

Some  characteristic  frequency. 


If  the  energy  in  each  eigenmode  is  approximately  equal,  the  following 
relationship  can  be  derived  from  the  above  composite  Q  equation: 


This  appears  to  be  the  approach  used  in  NBS  Technical  Note  1066  [Ref.  4] .  By 
considering  the  wave  number  -  w^/c  to  be  a  continuous  function  of  three 
continuous  variables,  m,  n,  and  p  which  correspond  to  the  integers  m,  n,  and 
1  in  Equation  4  above,  they  construct  a  lattice  sum  of  individual  Q  values 
for  a  wave  propagating  in  the  z  direction  within  a  cavity  (as  oriented  in 
Fig.  4).  This  is  accomplished  by  integrating  over  the  volume  of  a  spherical 

X  Z 


b 

,  Figure  4:  Cavity  orlencaclon  used  to  derive  fields 
(from  Ref.  4) . 
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shell  in  the  first  cctant  that  has  an  inner  radius  k  and  an  outer  radius 
k+Ak.  Thus  all  positive  values  of  m,  n,  and  p  are  included  that  satisfy  the 
equation: 


;  k  i  i  Jc+AJt 

N\  a  / 

1  b  I 

1  d  I 

(9) 

The  following  result  was  obtained  [Ref.  4]: 

0  .  2  - J - 

2  1.^(1  a,  11  (10) 

8Jc\a  D  dj 

[Note:  the  relative  permeability  of  the  cavity  walls,  Mr*  was  added  to 
correct  a  typo  and  the  variable  "c"  as  presented  in  the  reference  was  changed 
to  “d"  to  avoid  confusion  with  the  speed  of  light  constant  used  above.) 

where 

V  -  Volume  of  Che  cavity, 

S  -  Total  Internal  surface  area  of  the  cavity  walls. 

5s  -  Skin  depth  of  cavity  walls. 

Equation  10  appears  not  to  be  usable  for  the  definition  of  composite  Q 
desired  for  this  experiment.  This  is  primarily  because  only  two  sets  of 
modes  were  used  in  the  calculation  and  the  zero  order  modes  seem  to  have  been 
given  undue  weighting.  These  considerations  will  be  discussed  as  part  of  the 
following  analysis. 

To  completely  describe  the  fields  within  a  rectangular  cavity,  all  six  types 
of  modes  (TE‘*’,  TM<*>,  TE'*’,  TM'*>,  TE<y»,  and  TM<»>)  must  be  considered.  The, 
composite  Q  will  then  be  calculated  by  substituting  the  power  lost  to  the 
walls  and  the  energy  stored  by  each  mode  into  Equation.  7.  Then,  an 
expression  will  be  derived  for  the  net  Q  vrtilch  will  include  power  losses  to 
receiving  antenna  and  B-dot  probes.  With  minor  modification,  the  TE**’  and 
TH<*>  modes  from  NBS  Technical  Note  1066  [Ref.  4]  can  be  de.seribed  as 
presented  in  Equations  11  and  12. 
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TE‘*’  Modes: 


S,  =  0 

=  ky  COS  (k^)  sin  (kyy)  sin (k^z) 

Ey  =  -k^  sin  ik^)  cos  (kyy)  sin  (k^z) 

i/ {jr\kj)  kly  cos  (kyX)  cos  (kyy)  sin  {k^z) 

=  -1/  (jr\k^)  k^  iCj  sin  (k^)  cos  (kyy)  cos  (k^z) 
Hy  =  -1/  (jr\kj)  ky  JCj  cos  (k^)  sin(kyy)  cos  (k^z) 


TM<*’  Modes: 


^,  =  11/  (jkr)  kly  sin(k^)  sin  (kyy)  cos  (k^z) 

E^  =  -y\/ (jkj)  kj,  kg  cos  (kgX)  sin(kyy)  sin(KgZ) 

Ey  -  -i\/ (jkg)  ky  kg  sin  (kyX)  cos  (kyy)  sin  (kgZ) 

Hg^  0 

Hg  =>  ky  sin(kgX)  COS  (kyy)  COS  (kgZ) 

Hy  =  -JCj,  coslkgX)  sin(kyy)co3(kgZ)  (12) 


where  k,,  k,,  k^,  k,,  and  k^y  have  the  same  significance  as  in  Reference  4  and 
eta  is  the  intrinsic  impedance  of  the  medium  filling  the  cavity.  That  is: 

I-  iTfTt  nit  T-  iit 

kg  -  — ,  ky  »  -g-,  kg  -  -g-, 

/  kg  »  ^kg*ky*kg  (13) 


The  X  and  y  modes  can  be  easily  determined  by  noting  that  when  the  enclosure 
is  rotated  as  shown  ip  Figure  5,  the  following  transformations  occur: 


z  to  X  modes 
z-x 
x-y 
y-z 


z  to  y  modes 
z-y 
y-x 
x-r 


(14) 


Applying  Equation  14  to  Equations  11  and  12,  the  other  TM  and  TE  modes  are: 
TE‘*>  Modes: 


Eg*  Q 

Ey  *  kg  sin(kgX)  COS  (kyy)  sin(kgZ) 

Eg  •  •‘ky  sin(kgX)sin(kyy) cos  (kgZ) 

Hg*  1/ (jy\kg)  klg  sin  (kgX)  cos  (kyy)  cos  (kgZ) 

Hy  •  (jx\kg)  ky  kg  cos  (kgX)3in{kyy)  cos  (kgZ) 

Hg  •  •‘X/ (j^kg)  kg  kg  cos(kkx)cos(kyy)sin{kgZ)  (15) 
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Figure  5.  Three  different  directions  of  propagation. 

Modes: 


kyt  co3lk^)ain(kyy)3iTi{kgZ) 

^y  *  Uky)  ky  k^  3in(k„)co3(kj^)  aixiik^) 

E^^-t\/(jky)  ky  kj^  siTi(kj^)3inlkyy)co3(k^) 

«,  *  0 

Hy  ^  k,,  co3{k^)  3in(kyy)co3ikgZ) 

=  -ky  cos  (k^)  cos  (kyy)  sin  (Jc,z)  (16) 


TE<y>  Modes: 


Ey  •  0 

Eg"  kg  3in{kgX)3ia{kyy)coa(kgZ) 

Eg  "  -kg  cos  (kgX)  sin  (kyy)  sin 

Hy  -  1/  (jr\kg)  kgg  cos  (kgX)  sin  (kyy)  cos  <ic,z) 

Hg  "  -1/ (jr\kg)  kg  ky  co3(kgX)  cos  (kyy)  3in(kgZ) 

Hg  "  -1/  (jr\kg)  kg  ky  sin  (kgX)  cos  (kyy)  cos  (*^)  (17) 
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TM'y’  Modes: 

By  n/  (j^i)  kit  sin  ( Jcpr)  cos  (ir^y)  sin  ( Jr,2r) 

By*  -n/O'*,)  kt  ky  ain{ktt)s^Mkyy)coB\ktZ) 

Bx  =  ^jky)  kt  ky  cos  (k^x)  sin(JCyy)  sin(Jr,2) 

Hy  M  0 

H,  *  ky.co3(k,^)cos(kyy)ain(ktZ) 

Hy  •  -ky  ain(ktX)coa(kyy)coaiktZ)  (18) 


where 

ky,  -  ^jkl*kl ,  k„  »  {ki*ki 


(19) 


Now  that  at  all  six  types  of  irodes  have  been  defined,  it  can  be  noted  that 
the  TE  and  TM  are  orthogonal  in  time.  That  is: 


T  T 

j  F™  •  dr.  •  !  5'™  •  dc  “  t)  (20) 

.  0  0 

where  T  is  the  period.  The  fields  for  the  three  directions  are  not,  however, 
orthogonal  for  a  given  polarization.  Thus,  when  calculating  the  power  lost 
to  the  wall.'!  and  the  energy  stored  in  the  cavity,  the  TE  and  TM  modes  may  be 
treated  separattly,  but  the  fields  for  the  x,  y,  and  z  directions  must  be 
added  vectorial ly  within  the  TM  or  TE  mode  set.  Hence: 

TE<»-y  *’  Modes: 

B,  ^  (ky-k,)  cos(J(jcX)sin(Jlr^)8in(ir,z) 

By  •  (k,-k,)  sin  cos  (/rp^)  sin  (Jc^) 

•B,  ■  (k,-ky)  axn(k^)ain{kyy)caa(k,z) _ 

B,  •  {jy\kj)  {kl,~k,ky’‘kJc,Y.aimk,x)coa{kyy)  cos  ik,z) 

Hy^l/ijT\k,)  {klt-kyky-kyk,)  cos{k,x)ain{kyy)coa{k,z) 

H,»l/ijr]k,)  (kly~k,k,^kyk,)  coa{k,x)coa(kyy)ain{k,z)  (21) 


TM<*'y-*’  Modes: 

S,  •  tj/  ijk,)  {kl,-kj^y~kxk^  cos (Jc^)  sin (A:^)  sin(Jc,z) 
By  •  i\/ ijk,)  {kl,-k,ky-kyk,)  sin (it^)  cos  [kyy)  sin {k,z) 
B,  *  r\/{jk,)  {kly-kjz,~kyk,)  sin(>c^)  sin(l:^)co3  {Jk,z) 
H,  •  {ky-k^  sin()rjr’^)cos(A:^)cos(lr,2) 

Hy^  (k,-k,)  coB{k,x)Bin{kyy)coe{k,2) 

H,  •  (/c,-Jty)  COB  { k ,x)  COB  ikyy)  sin  {k,z) 


(22) 
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The  power  lost  to  the  cavity  walls,  for  small  losses,  can  be  calculated  using 
the  perturbation  method  [Ref.  2] : 


dS 

s 


(23) 


where  Rs  is  the  surface  resistance  of  the  cavity  walls  and  Ut  is  the 
components  of  H  tangent  to  the  wall  surfaces.  Applying  Equation  23  to 
Equations  21  and  22,  the  power  lost  at  all  six  walls  to  the  TE  and  TH  modes 
is : 

[  (ab+ad)  (ab+bd)  (k^-kyk^-kyk,) » 


■Pjw  • 


*  (ad^bd)  (kly-kjc,-kyk,)  >] 

[  (ab+ad)  (ky~k,)^*-(ah*bd)  (k^~k^)^*  (adi-bd)  (k^-ky)>] 


(24) 


kz 


ky 


Figure  6.  Surface  of  constant  k. 


The  different  amplitudes,  Afg  and  At^.  for  the  TE  and  TM  modes  were  added  to. 
indicate  that  the  asiplitudes  of  the  polarizations  as  defined  above  may  not  be 
equal.  Now  the  k,,  k^,  and  k,  dependency  can  be  reduced  to  strictly  a  k, 
dependency  by  assuming  that  the  lattice  spacing  is  small  compared  to  k^  and 
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vhere  ft,  and  5,  are  respectively  the  relative  permeability  and  skin  depth  of 
the  walls. 


Even  though  the  variation  of  k,  over  the  interval  k  S  k,  S  k+Ak  has  not  yet 
been  considered,  as  promised  earlier,  the  Justification  for  not  treating  the 
zero  order  modes  separately  (as  was  the  case  in  Ref.  4)  can  now  be  provided. 

Consider  in  the  equations  for  P*  (Eq.  24)  and  V  (Eq.  27),  the  integration 
over  the  surface  was  completed  in  spherical  coordinates  with  the  limits  of  8 
and  4  being  [0,ir/2].  The  zero  order  modes  lie  on  the  very  edges  of  the 
surface  of  Integration.  It  was  assumed  upfront  that  the  mode  density  was 
high  enough  to  consider  the  surface  continuous.  Therefore,  in  the  limit,  the 
zero  order  modes  contribute  nothing  to  the  surface  Integral.  In  Reference  4 
the  zero  order  modes,  even  though  there  arle  significantly  less  of  them,  are 
given  equal  weighting  in  the  sum  of  the  reciprocals.  To  be  less  biased,  each 
1/Q  value  should  have  been  multiplied  by  the  total  count  of  lattice  points 
contributing  to  that  1/Q  and  the  sum  then  divided  by  the  Count  of  all  lattice 
points  in  the  shell. 

Now  consider  the  case  when  the  source  pumping  the  cavity  has  its  power  spread 
out  over  some  spectrum.  It  is  convenient  to  express  this  type  of  signal  as  a 
power  spectral  density.  Let  the  source's  power  spectral  density  be  Pt(u)  in 
U/Hz.  Then  the  total  average  power  and  the  average  power  in  small  interval 
can  be  found  using  standard  communication  theory  [Ref.  5]: 
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^iC 

Au  <  u 


For  an  empty  cavity  (no  loads  or  receiving  devices) ,  the  power  absorbed  by 
the  walls  is  Identical  to  the  net  power  being  transmitted  into  the  cavity. 
Thus  the  power  lost  to  the  walls,  in  a  small  frequency  Interval  is  P^. 
Therefore,  the  energy  in  a  small  frequency  interval,  Au  can  be  found  by 
equating  Equations  29  and  30  and  substituting  into  Equation  6; 


Now  the  composite  Q  is  determined  by  svibstitutlng  the  above  energy 
expression  into  Equation  7  and  replacing  the  summations  with  integrals: 

CW  ■  - ^ 

du 

Clearly,  if  Px(u)  is  any  arbitrary  function,  then  the  composite  Q  will  be 
highly  dependent  upon  it.  But,  if  instead  the  spectrum  of  Px(w)  is  limited 
to  a  narrow  bandwidth,  BW  so  that  BW  «  thep  it  can  be  easily  shovm  that 
Q«qv  *  Qkr*  First  consider  a  definite  integral  over  a  small  interval  about  x, 
of  a  generic  function  f(x}  which  can  be  expressed  as  a  convergent  power 
series: 


j  £  (x)  dx  ■  j  dx 


jt,-*/*  I*'* 


♦  1  *  2x 
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The  bracketed  quantity  in  the  last  part  of  Equation  33  can  be  expanded  as  the 
difference  of  two  binoaial  series  with  all  the  even  terms  canceling.  Then, 
by  discarding  all  but  the  first  term  remaining,  the  equation  can  be 
simplified  to: 


Xc*t/2 

fix)  dx  = 


Applying  this  result  to  Equation  32,  the  following  result  occurs: 


J  U  ^p^lii)  d<J  =  B(f  Uc^PclO^) 

and 

j  PfMdct  •  BW  P^lu^) 

•  e-V/J 

hence 

.2 _ .0* 

”  4  S  Rg  2  5 


Xc-^/2 


1*0 


1  +  1 


Jf,. 


=  5  f  (Xg) 


(34) 


I 

I 


(35) 


The  error  introduced  by  discarding  all  but  the  first  term  of  the  binomial 
expansion  is  highly  dependent  on  the  even,  higher  order  derivatives  of  f(x), 
i.e. : 


Ezzor  * 


m 

E 


1 

2*-  (2ml)  I 


a  -  f  (x^)  »-* 


(36) 


Hence,  as  expected,  the  smoother  f(x)  Is  in  the  vicinity  of  x^,  the  better 
Equation  35  will  approximate  the  actual  integral. 

It  should  also  be  pointed  out  that  the  relative  permeability  for  most 
materials  (including  the  steels  commonly  used  to  construct  mode-stir 
ch^unbers),  except  ferrites,  is  »1  above  medium  frequencies  (300  KHz  to  3  MHz) 
[Ref.  6] .  Since  proper  mode  stirring  requires  a  high  mode  density,  useful, 
lower  frequency  bounds  range  in  the  hundreds  of  megahertz  for  a  typically 


16 


PL-TR-91-1036 


sized  chamber.  Therefore,  for  most  practical  purposes,  the  relative 
permeability  in  Equations  29  and  35  can  be  assumed  to  be  unity. 

Thus  far  only  empty,  vinloaded  cavities  have  been  considered.  Normally,  when 
the  Q  of  a  cavity  is  being  measured,  loads  in  the. cavity  should  not  be 
considered  separately.  However,  experience  has  shown  [Ref.  7]  that  at  the 
lower  end  of  the  spectrum,  theoretical  calculations  for  the  Q  always  tend  to 
drastically  overestimate  the  actual  Q.  A  large  part  of  the  loading  effect 
that  takes  place  can  be  attributed  to  the  antennas  used  to  either  pump  the 
cavity  or  to  measure  the  effective  power  density  within  the  cavity.  If  the 
actual  net  pumping  power  is  measured  (incident  minus  reflected)  and  used  as 
the  reference  for  the  Q  calculation,  then  the  transmit  antenna  can  be 
ignored.  However,  the  receiving  antenna  must  still  be  accounted  for. 


The  effect  of  a  receiving  antenna  (horn  antenna)  on  the  Q  of  a  cavity  can  be 
calculated  in  a  straightforward  manner.  First,  the  expected  value  of  the 
gain  of  the  antenna  in  a  mode-stir  chamber  is  1.  This  is  intuitively  obvious 
from  the  fact  that  all  the  energy  radiating  from  an  antenna  must  end  up  on  a 
closed  sphere  t.'.kt  surrovinds  it.  Additionally,  the  polarization  factor  of 
1/2  must  be  accounted  for  (Ref.  8|.  Therefore  the  effective  area  of  the 
antenna  is ; 


8x 


(37) 


The  effective  power  density  can  be  calculated  from  energy  density  in  the 
cavity  and  then  combined  with  Equation  3  to  determine  the  net  power  lost  to 
the  antenna: 

(38) 

2kf 


The  Q  of  a  cavity  with  perfectly  conducting  walls  but  loaded  with  an  antenna 
is: 


JiJL  . 
Plo.t 


(39) 
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The  same  approach  can  be  used  to  determine  the  effect  of  the  B-dot  probps  on 
the  overall  Q.  First,  however,  consider  that  a  B-Dot  probe  measures  only  the 
component  of  th®  B  field  normal  to  the  loop  surface.  Therefore,  to  determine 
an  equivalent  B„^  to  use  in  Equation  3  consider  a  wave  incident  on  the  B-Dot 
at  some  random  angle,  B,  relative  to  the  normal  vector  of  the  probe’s  loop, 
then  integrate  over  all  angles  to  find  the  average  value  of  B„^.  Noting  that 
the  projection  of  the  B  field  on  the  normal  to  the  loop  is  related  to  the 
cosine  of  the  angle  between  them: 


n 

(B*)  [b^  cos^e.  dd  -  Ab* 


Thus  from.  Equations  3,  6,  and  40: 


(40) 


0„ 


100  V 


thus 


(41) 


where  N,  and  Nj,  are  the  number  of  standard  antennae  and  the  number  of  B-dot 
probes  respectively.  It  should  be  noted  that  the  reciprocal  of  the  stun  of 
reciprocals  is  appropriate  because  the  energy  is  shared  by  all  the  loads  in 
the  cavity.  This  is  similar  to  the  Q  analysis  used  to  determine  the  effect 
of  a  lossy  dialectic  oh  a  cavity  (Ref.  2]. 


Figure  7  shows  the  theoretical  Qs  due  to  the  lossy  cavity  walls,  a  receiving 
antenna,  and  four  0.125  cm^  per  half  B-dot  probes  for  the  NT  enclosure.  Even 
with  four  B-dot  probes  (eight  half  probes),  the  associated  probe  Q  has  little 
impact  on  the  net  Q.  On  the  other  hand,  at  the  low  frequencies,  the 
receiving  antenna  can  have  a  significant  impact  on  the  net  Q.  In  fact,  for  a 
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le»9 

1««8 

la  +  7 
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0 

l*  +  5 

le*  3 

t««a 

l**i  '  1«*0  le*l 

Frequency  (GHz) 

Figure  7.  Effect  of  probes  and  antenna  on  the  net  Q. 

good  quality  chamber,  Che  net  Q  measured  at  Che  lover  frequencies  will  be 
almost  exclusively  due  to  Che  antenna.  As  can  be  seen  from  the  figure,  one 
can  define  a  minimum  critical  frequency,  f«,  where  the  Q  curve  due  to  the 
antenna  and  the  Q  curve  due  to  the  lossy  vails  cross: 


Theoretical  0  For  NT  Enclosure 


(42) 


(Note  that  was  not  included  in  the.  above  expression  due  to  the  earlier 
discussion  Oa.'  it  being  unity  at  the  typical  frequencies,  >3  MHz,  used  in 
mode-stir  chambers  and  for  Che  typical  metals  used.) 

Thus,  well  below  this  frequency  the  ancenna(s)  will  dominate  the  net  Q,  while 
well  above  this  frequency,  the  Q  for  the  lossy  walls  will  dominate. 
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Figure  8  shows  the  predicted  value  of  the  Q,  based  upon  Equations  39  and  41, 
for  the  2.74  X  3.05  X  4.57  m  steel  NIST  chamber  as  well  as  the  actual 
measured  Q.  A  complete  description  of  this  chamber  and  how  the  measurements 
were  made  appears  in  Reference  1.  The  theoretical  model  would  appear  to  fit 
the  shape  of  the  experimental  data  very  well  all  the  way  up  to  18  GHz.  There 
does,  however,  seem  to  be  a  factor  of  2  difference  between  the  prediction  and 
the  model.  This  difference  £  years  both  in  the  region  where  the  antenna 
dominates  the  Q  and  where  the  walls  dominate  the  Q.  If  the  difference 
appeared  only  in  the  latter  region,  it  could  be  explained  away  by  assuming 
extraneous  loads  in  the  chamber  such  as  dust,  test  stand  material,  etc. 
Because  the  factor  of  2  appears  in  both  regions,  though,  the  source  of  the 
error  must  be  common  to  both  Q  predictors.  The  only  parameters  that  Q,  and 
have  in  common  are  the  chamber  volume,  the  source  frequency,  and  of 
course,  the  experimental  data.  It  should  also  be  pointed  out  that  above 
18  GHz,  the  experimentally  measured  Q  begins  to  rapidly  exceed  the  predicted 
Q.  Neither  one  of  these  discrepancies  are  currently  explainable  and  furth.sr 
study  is  required. 
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FREQUENCY  STIRRING 


PREFERRED  METHOD 

The  conventional  method  of  mode-stirring  uses  a  metallic  paddle  wheel  to 
continuously  or  incrementally  change  the  boundary  conditions  within  the 
reverberation  chamber  while  maintaining  a  constant  pumping  frequency.  In 
this  manner,  field  homogeneity  is  achieved  by  averaging  the  contribution  of 
many  different  eigenfunctions  to  a  given  point  in  the  chamber.  In  the 
frequency  stir  method,  the  boundary  conditions  are  maintained  constant  while 
the  frequency  is  allowed,  to  vary  over  a  narrow  interval  about  some  center 
frequency.  Then  the  contribution  of  each  eigenfunction  to  a  given  spatial 
point  in  the  chamber,  in  Chat  narrow  band  of  frequency,  is  averaged  to 
provide  the  field  homogeneity.  The  advantage  of  the  frequency  stir  method 
over  the  mode  stir  is  that  by  spreading  the  power  over  a  band  of  frequency, 
Che  eigenfunctions  corresponding  Co  Chat  band  can  be  stimulated 
simultaneously.  Thus,  because  th>‘.  energy  stored  in  the  fields  in  an 
arbitrarily  small  volume  centered  at  any  given  spatial  point  is  the  average 
of  all  Che  energies  in  the  individual  elgenmodes,  field  homogeneity  is 
achieved  on  a  real  time  basis  (this  assumes  that  Che  mean,  of  the  energy 
density  contributed  by  the  elgenmodes,  in  the  narrow  band  of  frequency,  to  a 
given, point  in  the  chamber  is  spatially  stationary).  Real  time  homogoneity 
can  drastically  reduce  test  times.  Even  if  individual  eigenfunctions  are 
changed  by  adding  a  test  object  or  moving  an  existing  test  object  within  the 
chamber,  the  mean  over  Che  frequency  band  will  remain  constant  (within  the 
limits  of  the  sampling  variance). 

There  are  any  number  of  techniques  chat  can  be  used  to  provide  the  source 
frequency  agility  required.  At  minim;uB,  however,  Che  method  used  should  meet 
the  following  criteria; 

•  The  power  spectrum  should  be  flat,  across  Che  agility  bandwidth  to 
simplify  data  interpretation. 
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•  The  signal  should  be  ergodic  (or  at  least  loosely  time  stationary 
bver  the  maximum  averaging  period  of  the  test  article  and/or  the 
power  measuring  instrument)  so  that  real  time  homogeneity  can  be 
achieved. 

•  The  center  frequency  and  bandwidth  agility  of  the  source  should  be 
variable  over  a  wide  parameter  set  to  facilitate  test  flexibility. 

•  The  average  output  power  should  be  continuously  variable  or 
variable  in  3  dB  or  less  increments,  from  1  V  to  at  least  200  V 
to  allow  for  a  wide  range  of  cavity  quality  factors  and  test 
article  upset  levels. 

The  method  chosen  for  this  experiment  uses  band  limited  white  Gaussian  noise 
(WGN)  up-converted  by  performing  double  sideband,  suppressed  carrier 
modulation  with  an  RF  signal  from  a  synthesized  sweeper.  The  output  is  then 
amplified  with  a  TUT  amplifier  to  meet  the  higher  power  level  requirements. 
Center  frequencies  between  0.5  and  18  GHz  and  frequency  agile  bandwidths  of 
10,  20,  SO,  and  100  MHz  can  be  easily  realized.  In  addition,  the  power 
output  from  the  sweeper  can  be  varied  in  1-dB  increments  from  -120  dBm  to 
»  +15  dBm,  Therefore,  for  all  practical  purposes,  the  dynamic  range  of  the 
technique  Is  identical  to  the  dynamic  range  of  the  amplifier.  Figure  9 
demonstrates  the  theoretical  spectrum  of  the  signal  at  each  stage  of  the 
signal  processing.  The  synthesized  sweeper  is  used  to  generate  the 
monochromatic  signal  which  determines  the  center  frequency  of  the  output 
while  the  WGN  source  and  low  pass  filters  provide  the  frequency  agility.  One 
of  the  convenient  properties  of  WGN  Is  that  It  will  maintain  its  statistical 
properties  after  having  passed  through  a  linear,  invariant  system  such  as  a 
low  pass  filter  (Ref.  5). 

Figure  10  shows  the  actual  output  at  a  fixed  center  frequency  of  1.5  GHz, 
immediately  after  up— conversion,  of  the  technique  as  measured  by  a  spectrum 
analyzer. 
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LVJ 


Figure  9:  Frequency  stir  spectrum 
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Figure  10.  Output  of  the  up*converters  for  a  center  frequency 
of  1.5  GHz  and  all  foui:  agility  bandwldths.  . 
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FIELD  HOMOGENEITY  PREDICTIONS 

If  one  assumes  that  the  energy  density  at  an  arbitrary  point  within  the 
reverberation  chamber,  due  to  any  given  eigenmode,  is  a  random  variable  (this 
assumption  is  based  on  convenience  not  fact),  e^  with  definable  statistics, 
then  the  variation  of  the  mean  Ej  of  a  sample  of  n  eigenmodes  can  be 
described  using  simple  sample  theory  [Ref.  9]: 


where  <7,  and  o,  are  the  standard  deviation  of  the  random  variables  e  and  E 
respectively  (it  is  asstuned  that  the  sample  size  N  is  small  compared  to  the 
total  sample  space).  The  sample  size,  N  can  be  determined  directly  from  the 
mode  density  (Eq.  5)  as  function  of  the  agility  bandwidth,  BU  and  the  center 
frequency,  f^.  Hence; 


c* 

2  fj 

2xV  BW 

(44) 


Further,  the  standard  deviation  of  e  can  be  approximated  by  considering  the 
energy  density  due  to  an  Individual  eigenfunction  to  be  a  random  variable 
which  is  a  function  of  a  random  spatial  cdordlnate  (X,Y,Z). 


Let  X,  Y,  and  Z  be  Independent  random  variables,  each  uniformly  distributed 
on  [0,1] ,  i.e.'i 


fx(x)  ■  1;  Osxsl 
f,(y)  •  1;  Osysl 
fx(z)  >  li  Oizil 


(45) 


Their  joint  distribution  is: 


t{x,y.z)  -  f,(x)f^(y)fj(z>  ■  1 


(46) 
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Thus,  the  point  (X,Y,Z)  will  represent  a  landom  point  within  a  one  cubic  unit 
of  space  with  the  probability  of  selecting  any  point  in  that  volume  being 
equal.  Now  consider  a  small  volume  inside  a  rectangular  cavity  in  which  a 
single  eigenfunction  is  stimulated.  Further,  let  the  spatial  coordinates  be 
chosen  within  the  small  volume  so  that,  the  volume  encompasses  exactly  one 
spatial  period  of  the  eigenfunction  in  a  unit  cube.  A  normalized  fiuiction, 
g(X,Y,Z)  expressing  the  energy  den.sity  pattern  as  a  function  of  the  random 
variables  X,  Y,  and  Z  is: 


g{X,Y.Z)  -  8in^|-yAj  ain^|-|yj  ain^|-|-z| 


(47) 


The  expected  value  and  variance  of  g  can  be  easily  determined  [Ref.  10]; 


XXX 

•  Jffglx,y»z)  t(x,y,z)  dx  dy  dz  • 
0  0  0 


o*,  «  E[{g(Jf,  y,Z)-E(gU,  r.Z)]}*]  « 


(48) 


After  scaling  the  variance  with  the  amplitude  of  the  energy  density,  the 
standard  deviation  of  a  single,  randomly  sampled  eigenfxinctlbn  is  determined: 


•  w  a . 


8  W, 


opt  nr 

2nf  vy  8 


(49) 


Substitution  directly  into  Equation  44  yields; 


o  -  ^  19c* 

*  16/*  N 


(50) 
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It  should  be  noted  that  the  Q  in  the  above  equation  should  be  considered  a 
net  Q  as  defined  earlier.  However,  by  substitution  of  only  the  theoretical  Q 
from  Equation  35,  the  deviation's  dependence  on  the  frequency  and  cavity 
dimensions  can  be  determined  for  the  ideal  cavity; 

3  Pf  19  Po  o 
*  32  pj.  V  BUT 

[The  conductivity,  a,  under  the  square  root  sign  should  not  he  confused  with 
the  standard  deviation  symbol . ] 

At  first  glance,  it  would  seem  that  an  easy  way  to  Improve  the  homogeneity  of 
the  fields  in  a  given  cavity  would  be  to  simply  increase  the  surface  area,  S, 
by  corrugating  the  walls  or  maybe  increasing  the  volume.  This  would  in  fact 
work,  but  at  the  expense  of  the  maximun  field  levels  achievable  with  a  given 
Pf.  Clearly,  though,  increasing  the  agility  bandwidth,  BV,  or  the  frequency, 
f,  will  improve  the  field  homogeneity  without  sacrificing  field  amplitudes. 
Although  Equation  51  is  useful  for  investigating  the  energy  density 
homogeneity* S'  frequency  and  bandwidth  dependence,  it  is  desirable  to  formally 
define  a  homogeneity  parameter  that  is  Independent  of  the  power  pumped  into 
the  cavity  and  thus  the  amplitude  of  the  energy  density  within.  A  simple 
definition  (and  one  that  is  commonly  used),  yet  one  that  can  easily  be 
determined  experimentally,  is  the  maximum  average  energy  density  divided  by 
the  minimum  energy  density  (or  power  density)  measured  in  the  cayity, 
expressed  in  decibels,  le.: 


if  »  10  log-il^  .  10  log— 


It  is  difficult,  however,  to  use  statistics  to  predict  maxima. and  minima  of 
data  sets  when  the  density  function  is  not  knoim.  There  are  techniques  that 
can  be  employed  to  calculate  the  density  function  from  Equation  47  and 
Equation  52,  albeit  most  are  very  difficult.  Therefore,  the  preferred 
technique  is  that  of  Monte  Carlo  simulation.  Using  Equation  47,  N  random 
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Tfble  1.  Standard  deviation  of  field  ratios  for 
different  values  of  N. 

Monte  Carlo  Simulation  of  Field  Homogeneity 


r:^ 

a 

N 

a 

N  . 

a 

2 

11.35 

16 

2.54 

350 

0.515 

3 

8.47 

18 

2.35 

600 

0.380 

4 

6.58 

20 

2.23 

850 

0.323 

5 

5.61 

25 

2.03 

1100 

0.286 

6 

4.72 

29 

1.84 

1350 

0.260 

7 

4.49 

35 

1.65 

1600 

0.238 

8 

3.88 

45 

1.46 

1850 

0.220 

10 

3.47 

57 

1.26 

2100 

0.200 

12 

3.15 

82 

1.02 

2350 

0.190 

_ ilJ 

2.84 

100 

0.934 

2600 

0.186 

field  values  can  be  produced  and  averaged  by  using  three  randomly  generated 
values  for  X,  Y,  and  Z.  This  is  done  twice,  and  the  ratio  is  saved  in 
decibels.  This  is  repeated  for  a  large  number  of  points  (2000'-5000  in  this 
simulation)  and  then  the  standard  deviation  is  calculated  using  the  usual 
methods..  Table  1  lists  the  standard  deviation  for  several  different  values 
of  N.  Not  knowing  the  density  function  for  Equation  52,  relating  the 
standard  deviation  in  Table  1  to  Equation  52  is  not  possible.  If  it  is  ' 
assumed  a  priori  that  the  data  are  normally  distributed  (evidence  of  this 
will  be  provided  later) ,  then  a  rule  of  thumb  criterion  can  be  established. 
For  example,,  to  determine  a  level  in  which  99  percent  of  the  magnitude  of 
the  ratio  data  will  be  >  1^0.99.  2  1/2  standard  deviations  can  be  used 
(»  99  percent  of  the  area  under  the  normal  distribution  falls  between 
±  2  1/2  standard  deviations)  to  predict  the  experimental  results.  The 
selection  of  the  99  percent  threshold  is  arbitrary,  but,  as  will  be  seen  in 
the  results,  seems  to  relate  the  measured  standard  deviation  to  the  measured 
homogeneity  very  well,  l.e.,  ^^•2.5  a.  For  example,  if  one  desires  a 
homogeneity  of  3  dB,  then  from  Table  X,  at  least  57  modes  would  need  to 
stimulated  simultaneously. 
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Now  the  agility  bandwidths  required  to  ^»chieve  the  same  field  uniformity 
results  as  achieved  by  the  conventional  paddle  wheel  technique  used  in  a  full 
size  chamber  such  as  the  NIST  chamber  mentioned  in  the  introduction  can  be 
predicted.  As  stated  earlier,  field  uniformities  of  <  +/-  3  dB  above  1  GHz 
and  <  +/-  2  dB  above  2  GHz  were  achieved  for  the  40  m^  chamber.  Using  the 
99  percent  rule  of  thumb  discussed  above,  57  modes  would  be  required  for 
above  1  GHz  and  somewhere  in  the  neighborhood  of  200  modes  would  be  required 
above  2  GHz.  Applying  Equation  3,  the  agility  bandwidths  can  be  calculated 
as  s  1.5  MHz  and  1.3  MHz  for  the  1  GHz  and  2  GHz  range,  respectively.  Of 
course  this  is  an  ideal  situation,  and  therefore  a  more  reasonable  estimate 
is  probably  in  the  area  of  5  MHz . 


COUPLING  CROSS  SECTIONS 


One  of  the  fundamental  motivations  for  conducting  low  power  microwave  tests 
is  to  determine  inexpensively  a  system  or  subsystem's  coupling  cross  section 
or  effective  area.  Cross  section  is  defined  as  [Ref.  11]; 


a(f) 


(53) 


where ; 


<T(f)  Coupling  cross  section. 


P{,  >  Power  delivered  to  the  measuring  instrument. 

Hv(f)  «  Voltage  transfer  function  (voltage  impressed  at  the  test 
point  divided  by  the  magnitude  of  the  incident  E  field) . 


^  Input  impedance  of  the  measuring  instrument. 


5iiic(f)  “  Power  density  incident  on  the  system. 

The  first  observation  that  can  be  made  about  Equation  53  is  that  it  assumes 
that  the  test  object  is  being  radiated  with' a  plane  wave.  Secondly,  it  also 
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assumes  that  the  source  of  radiation  Is  at  a  single  frequency.  Neither  of 
thesf-  assumptions  are  correct  in  the  case  of  frequency  stirring.  As 
discussed  earlier,  a  B-dot  probe. inside  a  cavity  with  an  isotropic  field, 
will  measure  1/2  of  the  power  density  associated  with  a  wave  traveling  in  a 
single  direction  that  is  properly  oriented  to  the  probe.  If  you  consider  a 
point  inside  this  cavity  with  isotropic  fields,  then  the  power  density 
Impinging  on  that  point  will  continuously  change  direction,  yet  the  average 
magnitude  of  the  power  density  will  remain  constant  and,  as  mentioned,  will 
be  twice  that  measured  by  a  B-dot.  For  lack  of  a  better  term,  this  power' 
density  wij.1  be  called  fringe  power  density  (borrowed  from  microwave  oven 
terminology).  Now  consider  the  surface  of  a  large  test  object;  the  power 
incident  on  the  surface  is  now  limited  to  2it  steradians.  If  one  is  simply 
interested  (when  studying  heating  effects,  for  example)  in  the  average  power 
density  falling  on  the  surface  (call  it  surface  power  density)  it  can  be 
calculated  as: 


%/i 

P.  •  (  P*  cosO  d6  ■  —  P,  ^54) 

0 

Where  P«  and  P{  are  the  surface  and  fringe  power  densities  respectively. 
However,  when  studying  the  coupling  due  to  a  port  of  entry  (POE),  the 
expected  value  of  the  directivity  will  always  be  one.  However,  the 
polarization  factor  will  vary  from  POE  to  POE. 

In  the  case  of  subsystem  which  would  norjMlly  exist  In  an  equipment  bey  or 
cavity,  it  makes  sense  to  use  the  fringe  power  density  when  determining  cross 
sections  because  this  is  typically  what  is  used  when  measuring  the  shielding 
effectiveness  of  the  full-up  system.  In  other  words,  the  full-up  object  is 
illuminated  with  a  plane  wave  while  the  fringe  power  density  is  measured 
within  the  cavity.  The  ratio  of  the  two  then  relates  the  fringe  power 
density  within  the  cavity  to  the  plane  wave.  Thus,  cross  sections  measured 
relative  to  the  fringe  power  density  can  be  related  directly  to  power  density 
incident  on  the  full-up  system.  A  question  immediately  arises,  however,  if 
the  coupling  cross  section  at  a  given  circuit  node  within  a  full-up  system  is 
measured  in  a  reverberation  chamber,  how  can  this  information  be  related  to 
plane  wave  illumination  in  an  anechoic  chamber,  for  example?  If  the  cross 
section  for  a  full-up  system  is  measured  relative  to  the  fringe  power 
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density,  then  this  is  equivalent  to  radiating  the  test  object  in  an  anechoic 
chamber  with  a  power  density  at  target  equivalent  to  the  fringe  power  density 
from  many  different  directions  with  many  different  polarizations  and  then  ■ 
averaging  the  results.  Another  useful  way  of  interpreting  the  cross  section 
as  measured  in  chamber  wo\'.ld  be  as  the  expected  value  of  the  cross  section 
that  would  occur  if  the  target  was  shot  at  from  a  random  direction  with  a 
random  polarization. 

So  far  only  the  effect  of  isotropic  illumination  on  the  cross  section  has 
been  discussed.  The  effect  of  spreading  the  spectrum  of  the  pumping  source 
will  now  be  considered. 

Because  the  measuring  instrument  effectively  measures  the  total  average  power 
being  received  by  the  probe  over  the  agility  bandwidth,  measuring  the  cross 
section  using  this  method  is  equivalent  to  cbnvoluting  a  gating  function  with 
the  narrow  band  cross  section  and  then  dividing  by  the  fringe  power  density, 
i.e. : 


(55) 

This  has  the  same  effect  as  mathematically  smoothing  <7(f)  over  the  same 
bandwidth.  The  primary  difference  is,  however,  that  the  smoothing  is  done  ,in 
a  real  time  manner  and  is  independent  of  the  sampling  resolution.  In  other 
words,  if  the  sample  interval  was  greater  than  the  smoothing  window, 
mathematical  smoothing  would  not  Votk  and  this  method  would. 

As.  shown  in  Figurefs  11  and  12,  details  in  the  crosk  section  that  span 
frequency  intervals  less  than  the  agility  bandwidth  will  tend  to  be  washed 
out,  while  those  that  span  greater  intervals  will  remain.  In  fact,  when  the 
bandwidth  of  the  resonance  is  less  th/^n  the  bandwidth  of  the  source  pumping 
the  cavity,  the  resonance  will  be  reduced  by  the  ratio  of  the  resonance  to 
the  source  bandwidth.  This  may  create  a  problem  if  the  test  asset  has 
coupling  resonances  with  bandwldths  smallbr  than  the  agility  bandwidth 
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Figure  11.  Measured  cross  section  when  the  bandwidth  of  resonance 
Is  smaller  than  the  bandwidth  of  the  gating  function. 
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Figure  12.  Measured  cross  section  when  the  bandwidth  if  resonance 
is  greater  than  the  bandwidth  of  the  gating  function. 
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required  to  maintain  adeqtiate  stirring.  Typically,  however,  test  object 
resonances  less  than  a  few  tens  of  megahertz  are  due  to  the  modes  set  up  in 
the  test  article  itself.  These  resonances  are  normally  of  little  interest 
because  of  their  dependence  on  the  finer  details  of  the  test  asset.  For 
example,  if  a  cable  run  on  one  aircraft  is  slightly  different  than  another  of 
the  same  type,  drastic  differences  can  b^  realized  in  the  details  of  the 
cross  section  even  though  the  trend  data  remain  the  same.  Objects  as 
insignificant  as  an  aluminum  soda  container  can  make  significant  differences 
in  the  frequencies  at  which  these  narrow  resonances  occur.*  Except  for  the 
introduction  of  some  ambiguity  at  the  band  edges,  the  wider,  POE  resonances 
will  pas-s  through  this  technique  unscathed,  as  can  be  seen  from  Che  amplitude 
of  Che  resulting  convolution  in  Figure  12. 

SHIELDING  EFFECTIVENESS 

Shielding  effectiveness  (SE)  is  a  dimensionless  quantity  which  relates  the 
fields  measured  inside  a  test  object  to  chose  that  would  be  measured  if  the 
test  object  were  not  present  and  can  be  experimentally  measured  by  simply 
measuring  the  fields  in  the  test  object  using  a  probe  and  Chen  dividing  Chat 
by  the  free  field  measurement  made  with  Che  same  probe.  VThen  doing  this 
measurement  using  Che  frequency  stir  technique,  the  considerations  that  were 
made  for  Che  coupling  cross  section  measurement  must  also  be  made  for  Che 
shielding  effectiveness  measurement.  That  is,  Che  resulting  information 
obtained  is  an  average  or  expected  value  over  many  angles  of  incidence  and 
polarizations,  and  the  data  are  convoluted  with  Che  pumping  source  function.. 

UPSET  THRESHOLD  TESTING 

One  of  the  most  useful  applications  of  reverberation  chamber  testing  is  its 
use  to  determine  the  threshold  for  electromagnetic  upset  of  a  unit  under 
test.  This  is  primarily  because  very  large  field  levels  can  be  achieved  in 
these  chambers  using  modestly  low  power  driving  sources.  In  the  conventional 
paddle  wheel  mode-stir  chamber,  upset  testing  can  be  a  lengthy  process 


♦Price,  R.,  ,ec  al.  The  Pepsi  Syndrome  as  it  Applies  to  High  Power  Microwave 
Susceptibility  Analysis.  JAYCOR.  Albuquerque,  NM,  August  1990  (Unpublished) . 
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becjiase  for  each  power  level  attempted,  the  paddle  wheel  must  be  rotated 
through  many  positions  at  a  rate  less  than  the  test  object's  response  rate. 
Once  the  lowest  powSr  level  in  which  an  upset  occurs  is  recorded,  the  power 
density  fn  the  chamber  can  be  measured.  The  difficulty,  however,  in 
interpreting  these  results  is  that  when  the  upset. occurs,  there  is  at  no  time 
a  truly  uniform  field  within  the  chamber.  In  other  words,  the  front  of  the 
test  object  may  be  in  a  hot  spot  and  the  rear  of  the  object  may  be  in  a  cold 
spot  when  the  upset  takes  place.  The  frequency  stir  method  eliminates  this 
ambiguity  by  providing  a  real  time  homogenous,  isotropic  field.  Also, 
because  the  field  is  the  same  throughout  the  chamber,  the  power  density  can 
be  monitored  while  the  source  power  is  slowly  raised  until  an  vtpset  occurs. 

There  is  cne  point  that  should  be  noted  when  using  the  frequency  stir 
technique  to  cause  upset.  That  is  the  fact  that  the  signal  being  received  by 
the  imit  under  test  contains  up— converted  noise  that  when  detected  by  the 
test  object  may  be  in-band  of  the  object's  electronics.  In  the  paddle  wheel 
method,  the  detected  signal  ‘s  Just  a  DC  level.  The  unit  under  test  may 
respond  differently  to  a  noisy  signal  than  to  a  DC  signal.  The  frequency 
stir  upset  method  was  successfully  used  on  an  aircraft  subsystem.  Testing 
was  completed  in  <5  minutes.  The  results  of  the  test  are  not  Included  here 
due  to  their  security  classification  and  the  fact  that  there  currently  are  no 
other  data  available  on  this  device  for  comparison. 
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TEST  CONFIGURATIONS 

FIELD  UNIFORMITY  SWEPT  MEASUREMENTS  YNT  Enclosure) 

Figure  13  shows  the  general  test  setup  used  to  conduct  the  field  uniformity 
tests  for  the  NT  enclosure.  Four  B-dot  (JAYCOR  TLB- 3,  8  GHz,  0.125  cm^/half) 
probes  were  placed  at  different  locations  within  the  reverberation  chamber. 

As  depicted  in  Figure  14,  they  were  placed  in  random  locations  with  all, three 
orientations  (x,y,  and  z)  represented.  A  three-dimensional  schematic  of  the 
NT  reverberation  chamber  appears  in  Figure  1.  As  shown,  it  is  a  6  x  4  x  5  ft 
welded  aluminum  enclosure  with  three  2  x  2  ft  access  panels  with  RF  gaskets 
and  secured  with  bolts  every  3  in  around  their  perimeter.  The  microwave 
launch  into  the  cavity  is  accomplished  using  a  broad  band  (1  to  18  GHz)  horn 
positioned  at  nonorthogohal  angles  to  the  cavity  wall  to  facilitate  multiple 
mode  excitation.  Equipment  control  and  data  acquisition  are  facilitated  via 
a  486  computer  using  a  National  Instruments  GPTB  board.  The  HP83620A 
Synthesized  Sweeper  provides  the  user  selected  frequency  range  and  power 
level  to  the  up-converter  bank.  Though  not  shown  in  Figure  13,  external 
level  control  was  provided  via  a  coupler  and  HP423B  Crystal  Detector  placed 
at  the  output  of  the  TVT  bank.  External  leveling  was  used  to  maintain  the 
TUTs  at  an  output  power  level  of  »43  dBm.  This  helped  to  minimize  TUT 
harmonic  distortion  and  protect  other  equipment  in  the  system  from  being 
damaged  by  excessive  power. 

The  frequency  agility  is  provided  by  the  NC7907  Noise  Source.  The  noise 
source's  bandwidth  and  output  power  can  be  set  by  Che  486  computer.  In  order 
to  vary  the  maximum  100-MHz  bandwidth  output,  the  noise  source  has  four  user 
selectable  internal  low  pass  filters  (5,  10,  25,  and  50  MHz).  Further,  the 
maximum  output  power  of  30  dBm  can  be  varied  via  an  internal  programmable 
attenuator.  The  control  software  developed  for  this  experiment  automatically 
selects  Che  appropriate  filter  and  attenuator  value  based  upon  the  desired 
agility  bandwidth  for  the  test  at  hand.  The  up-converter  bank  provided 
double  sideband,  suppressed  carrier  modulation  of  the  microwave  carrier  by 
the  band  limited  WGN.  Due  Co  the  double  sideband  nature  of  the  modulation, 
the  output  of  the  up -converter  bank  has  twice  the  bandwidth  of  the  output  of 
Che  noise  source. 
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Figure  13.  Experimental  setup  for  swept  tests  on  the  NT  enclosure. 
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Figure  14.  The  NT  enclosure:  view  of  horn  and  prohes. 
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Output  from  the  ,B-dot  probes  is  converted  to  a  signal  measurable  by  the 
HP8757C  Scaler  Network  Analyzer  via  one  of  four  HP85025A  AC/DC  Detectors 
(0.01  to  18  Ghz,  23  dBm  max).  The  HP8757C  allows  for  a  great  deal  of 
flexibility  in  making  measurements  of  incoherent  signals  while  maintaining  a 
high  degree  of  measurement  accuracy.  Simple  data  analysis  such  as  detector 
ratios  and  averaging  can  be  performed  in  near  real  time  during  the  actual 
measurement  by  the  scaler  analyzer.  The  computer  controller  breaks  the  user 
supplied  start  and  stop  frequencies  into  individual  bands  that  match  the 
bands  of  the  TUTs  and  up -converters,  and  then  it  instructs  the  HP8757C  to 
acquire  the  appropriate  data  for  each  band  while  it  provides  the  switching 
signals  to  the  TWT  and  up-converter  banks.  The  data  from  each  band  are  then 
sampled  in  accordance  with  the  user  directed  nvunber  of  data  points  and  a 
complete  sweep  is  constructed. 

FIELD  UNIFOPwMITY  SWEPT  MEASUREMENTS  (GIFT  Box) 

The  GIFT  Box  is  a  30  x  36  x  36  in  welded  aluminum  enclosure  with  one  large 
(27  in  diameter)  access  hatch  in  the  front  and  two  other  smaller  (14.5  in 
diameter)  access  panels,  one  on  the  side  and  one  on  the  back.  In  addition, 
as  shown  in  Figure  2,  there  is  one  4.5  in  diameter  feed-through  panel  on  the 
front  where  the  signals  for  the  two  B-dot  probes  were  routed  through. 

Figures  15  and  16  depict  the  test  configuration  used  during  the  GIFT  Box 
swept  field  uniformity  tests.  With  the  exception  of  the  amplifiers  and  the 
number  of  probes,  this  test  configuration  is  very  similar  to  that  used  for 
the  NT  enclosure  tests.  The  TVTs  were. not  available  at  the  time  this  test 
was  conducted;  instead,  a  set  of  high  gain  small  signal  amplifiers  was  used 
on  the  receiving  side  of  the  circuit  to  boost  the  B-dot  output  to  levels  well 
within  the  dynamic  range  of  the  Scaler  Analyzer. 

IHS  ■a.J1EAgVRgMOTg.jCMI..3Ti<l,C.III.Joxl 

The  test  configuration  used  to  measure  the  Q  of  the  NT  enclosure  is  shown  in 
Figure  17  (the  GIFT  box  setup  is  similar  except  small  signal  amplifiers  were 
used  rather  than  TWTs).  The  significant  difference  between  the  two 
configurations  is  the  same  as  in  the  case  of  Che  swept  tests  (TVT  amplifiers 
versus  small  signal).  Further,  computer  control  of  the  sources,  modulators, 
and  amplifiers  is  similar  to  the  swept  ease  except  that  a  single  frequency  is 
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sele<;ted  for  each  measurement  and  the  S3mthesized  Sweeper  is  placed  in  the 
pulse  modulation  mode.  The  HP3314A  Function  Generator  is  used  to  supply  a 
TTL  compatible  square  wave  to  the  Synthesized  Sweeper  to  chop  the  output  of 
the  sweeper  at  a  relatively  low  frequency  (10  to  100  KHz) .  The  chopped 
signal  is  then  supplied  to  the  up -converter  bank  for  modulation  by  the  band 
limited  WGN  source.  Both  the  horn  input  signal  and  the  signal  received  by 
the  B-dot  probe  are  detected  using  HP423B  Crystal  Detectors  and  sent  to  the 
DSA  602  Digitizing  Signal  Analyzer  for  averaging  and  rise  time  measurements. 
It  should  be  noted,  that  because  negative  detectors  are  used,  the  B-dot  signal 
rise  time  is  equivalent  to  the  relaxation  time  of  the  enclosure. 


Figure  15.  Experimental  test  setup  for  swept  tests  on  the 
GIFT  box  enclosure. 
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Figure  16.  The  GIFT  box  enclosure:  view  of  horn  and  probes. 


Figure  17.  Test  configuration  used  to  measure  Q  of  the  NT  enclosure 
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TEST  PROCEDURES 


Field  uniformity  tests  were  conducted  on  both  the  NT  and  GIFT  Box  enclosures 
to  determine  the  effectiveness  of  the  frequency  stir  method  at  reducing 
localized  average  field  gradients.  Using  the  test  configurations  described 
above,  B-dot  probes. were  used  to  measure  the  power  ratios  of  the  magnetic 
fliix  densities  at  randomly  selected  locations  and  vector  orientations  within 
the  enclosures.  Using  Equation  3,  for  the  monochromatic  case  (no  modulation) 
and  with  the  frequency  fixed,  the  ratio  of  the  power  measured  at  the  analyzer 
between  any  two  probe  locations,  A  and  B,  can  be  expressed  as; 


^  Vi/(2z)  ^  ( Jodf  =  fisdf 

Pb  Vl/{2z)  l^oaj  l^oaj 


where 


(56) 


P^Pg  •  The  power  ratio  of  the  probe  outputs  as  directly  measured  by 
the  analyzer. 

Vg  -  The  amplitude  of  the  voltage  produced  by  the  B*dot  probes. 

z  Impedance  of  the 

Hoa/Hqb  “  The  ratio  of  the 

along  the  axis  of  tjhe  probes. 

When  modulation  is  used  the  signal  being  measured  by  the  analyzer  is  no 
longer  monochromatic  and  therefore  the  power  ratio  analysis  is  slightly  more 
complex.  Using  the  probe’s  transfer  function.  Equation  2,  the  power  measured 
by  the  analyzer  will  be; 


measuring  device.  . 

magnetic  fields  at  locations  A  and  B  oriented 
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p  =  1^  S(<i))  dta  (57) 

2-K  lit  ' 


where 

{i>e  -  Center  frequency  of  the  received  signal. 

BV  -  Bandwidth  of  the  received  signal. 

S(c>>)  ••  Power  Spectral  density  of  the  magnetic  field  at  the  probe 
location. 


Therefore,  the  ratio  of  the  power  measured  can  be  expressed  as: 

P  I* 

Pa  I*  S,(«)  flto 

Jit  t-aa/2 


(58) 


Clearly  this  equation  can  only  represent  the  power  ratio  of  the  magnetic 
fields  when  the  square  of  the  magnitude  of  the  transfer  function  C  Is 
constant  over  the  Interval  Ug-SV/2  <  la  <  ta^+iM/2.  From  Equation  2,  the 
percent  variation  of  can  be  easily  written: 


y 

Bztot 


iG(iw)  I*  ■  aI  »* 
Ay  ^  2  aI  <a  BW 

y  *  Ala* 


BW 

nf 


(59) 


At  the  lowest  frequency  (1  GHz)  and  the  widest  agility  bandwidth  (100  MHz) 
used  In  this  experiment,  the  variation  Is  *3  percent.  It  is  therefore 
evident  that  for  any  given  frequency  thb  field  variation  between  two  points 
can  be  determined  directly  from  the  test  configuration  discussed  earlier.  It 
Is  not,  however,  practical  to  fix  the  frequency  and  move  the  probes  around 
the  volume  of  the  enclosures  until  the  maximum  field  variations  are  found. 
Instead,  the  maximum  and  minimum. field  locations  are  forced  to  move  within 
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the  volume  by  sweeping  the  center  frequency  of  the  injected  signal.  As  the 
center  frequency  is  swept,  new  modal  frequencies  are  encompassed  in  the 
signal's  bandwidth  while  lower  medal  frequencies  drop  off.  Thus  the  maxima 
and  minima  of  the  fields  move  around  the  enclosure  passing  by  the  probes  as 
they  move.  The  spatial  variation  of  the  fields  can  therefore  be  approximated 
by  observing  the  maximum  variations  in  the  swept  data  provided  a  sufficiently 
large  frequency  interval  is  used. 

Using  the  method  described  above,  swept  tests  were  conducted  on  the  NT 
enclosure  from  1  through  4  GHz  with  both  monochromatic  and  up-converted  noise 
with  bandwidths  of  10,  20,  50,  and  100  MHz.  Ratio  measurements  were  taken 
directly  usivig  the  HP85025A  detectors  as  described  in  the  test  configuration 
above.  The  copper  jacket  0.141  cable  used  between  the  probes  and  the  feed¬ 
throughs  at  the  bulkhead  were  cut  to  the  same  length  and  hence  their 
insertion  losses  were  measured  to  be  within  a  few  tenths  of  a  decibel  of  each 
other  and  therefore  ignored.  The  GIFT  Box  sweeps  were  conducted  from  1 
through  8  GHz.  As  with  the  NT  enclosure,  both  monochromatic  and  up-converted 
noise  with  bandwidths  of  10,  20,  50,  and  100  MHz  were  used.  However,  because 
the  small  signal  amplifiers  wer6  used  on  the  receiving  side  of  the  circuit, 
separate  sweeps  for  each  probe  were  taken.  The  sofMare  was  then  used  to 
convert  the  measurements  to  a  ratio  similar  to  that  of  the  NT  enclosu.re  data. 


The  method  used  here  to  measure  the  Q  of  the  two  enclosures  is  a  modification 
of  the  method  developed  by  the  Naval  Surface  Weapon  Center.  The  I/e  decay 
time  of  the  energy  stored  in  the  cavity  can  be  expressed  as  [Ref.  12]: 


(60) 


The  Weapon  Center  measured  the  energy  decay  waveform  produced  by  a  pulsed  RF 
input  for  many  positions  of  the  paddle  wheel  of  a  mode-stir  chamber.  Then  by 
calculating  the  fall  time  of  the  ensemble  average,  they  determined  the 
composite  Q  from  the  above  equation. 
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Since  there  Is  no  paddle  wheel  vised  in  this  technique,  the  Q  of  the  cavities 
was  measured  slightly  differently.  Employing  the  frequency  stir  method  and 
the  test  configuration  described  earlier,  the  relaxation  time  of  the 
enclosures  was  determined  by  using  the  measured  fall  time  of  the  ensemble 
average  of  the  detected  pulse  train.  This  would  seem  to  be  a  valid  method  to 
determine  the  fall  time  of  the  stored  energy  in  the  cavities  because  the 
ensemble  defined  by  each  pulse  contains  information  from  both  the  noise 
envelope  and  the  energy  decay  waveform,  but  since  t  'e  ensemble  average  of  the 
noise  envelope  is  a  constant,  the  ensemble  average  of  many  of  the  pulses  is 
proportional  to  the  stored  energy  waveform.  Therefore,  the  fall  time  of  the 
ensemble  average  is  equal  to  the  fall  time  of  the  stored  energy  waveform. 
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SWEPT  MEASUREMENTS 

Figure  18  depicts  the  results  of  the  swept  measurements  for  the  NT  enclosure 
from  1  to  4  GHz.  As  Indicated,  agility  bandwldths- of  0,  10,  SO,  and  100  MHz 
are  shown  (the  results  of  the  20  MHz  agility  bandwidth  sweeps  are  not  shown 
because  of  their  similarity  to  the  10  MHz  data)  for  the  measured  power  ratio 
of  probes  C  and  D.  These  data  are  typical  of  the  data  taken  for  all 
combinations  of  the  probe  A.  B,  C,  and  D.  The  theoretical  number  of 
elgenmodes  stimulated,  predicted  standard  deviation,  measured  standard 
deviation,  and  ^  as  defined  by  Equation  52  are  given  In  Tables  2  and  3  for 
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Figure  18.  Measured  power  ratios  of  probes  C  and  D  for  the 
NT  enclosure  with  an  agility  bandwidth  of  0,  10, 
50  and  100  MHz. 


44 


PL-TR-91-1036 


1  GHz  frequency  ranges  for  each  agility  bandwidth.  The  Monte  Carlo 
simulation  predicts  the  standard  deviation  well  for  values  >1.  It  is 
believed  that  the  setup's  random  and  systematic  errors  begin  to  dominate  the 
standard  deviation  measurements  when  the  standard  deviation  begins  to  drop 
below  1.  Therefore,  data  above  4  GHz  are  not  shown  for  the  GIFT  box.  The 
largest  source  of  these  errors  is  probably  due  to  jthe  systematic  errors 
between  the  two  probes  used  for  the  measurement.  The  particular  probes  used 
here  have  shown  differences  as  large  as  a  few  decibels  as  the  frequency  and 
E-field  orientation  are  varied. 

As  mentioned  earlier,  a  normal  distribution  seems  to  be  an  adequate  modal  for 
the  density  distribution  of  the  ratio  data.  Figure  19  shows  the  typical 
results  of  fitting  the  ratio  data  (in  decibel  format)  to  a  normal 
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distribution.  The  histograms  represent  the  relative  frequency  distributions 
(RFD)  and  the  cumulative  relative  frequency  distributions  (CRFD)  for  the 
actual  data.  These  RFDs  are  formed  by  dividing  the  abscissa  of  the  data  into 
equally  sized  cells  and  counting  the  number  of  points  falling  in  each  cell. 
The  cell  count  is  then  divided  by  the  total  niomber  of  points.  The  CRDF  is 
simply  a  nmning  sum  of  the  cell  counts .  The  continuous  plots  are  the  RDF 
and  CRDF  of  the  theoretical  model  being  fitted. 

The  goodness  of  fit  test  used  is  the  Person's  chi-square  statistic  [Ref.  ’9]. 
Tbe  smaller  the  fit  value,  the  better  the  fit.  In  fabt,  by  looking  up  the 
fit  value  in  a  table  of  chi-square  distributions  of  n  degrees  of  freedom,  the 
probability  that  discrepancies  of  this  order  (the  fit  value)  or  larger  would 
be  seen,  given  the  model  is  correct,  can  be  found.  In  all  cases,  the 
distribution  of  the  ratio  data  (in  decibels)  proved  to  have  acceptable  normal 
distribution  fits. 


Figure  20  shows  the  quality  factor  of  both  enclosures  as  measured  using  the 
technique  described  earlier.  No  attempt  is  made  at  predicting  the  Q  based 
on  the  earlier  theoretical  Q  discussions  because  the  wall  conductivity  is 
difficult  to  assess  for  these  enclosures.  They  both  are  in  very  poor 
condition,  and  have  accumulated  scale  on  their  interiors. 


CONCLUSION 


From  the  data  gathered  to  date,  it  would  seem  that  the  frequency  stir 
technique  has  potential  application  in  electromagnetic  interference  and 
susceptibility  testing.  For  practically- sized  chambers,  field  uniformities 
equivalent  to  the  p&ddle  wheel  mode-stir  technique  can  easily  be  achieved 
with  relatively  narrow  agility  bandwidths  (5  MHz)  while  test  time  is 
dramatically  improved.  Another  potential  application  of  the  frequency  stir 
technique  is  to  directly  pump  a  large  test  article  with  the  up -converted 
noise  so  that  its  own  cavities  become  the  mode-stir  chamber.  This  may  be 
helpful  in  both  upset  testing  for  assets  too  large  to  fit  in  a  chamber  and 
sniff  testing  (finding  the  leaky  apertures  by  using  a  probe  to  sniff  around 
the  exterior  of  the  test  article).  Lastly,  because  of  the  inherent  smoothing 
of  the  coupling  dat,a  that  takes  place  when  using  this  technique,  it  may  have 
application  in  plane  wave  illumination  in  an  anechoic  chamber  to  acquire 
coupling  trend  data  directly  or  reduce  the  dynamic  range  required  by  the 
measuring  instruments. 

It  has  also  been  demonstrated  that  when  predicting  the  Q  of  a  reverberation 
chamber,  the  effect  of  the  measuring,  and  in  some  cases  the  transmitting 
antenna,  must  be  considered.  They,  in  effect,  act  as  apertures  in  the 
chamber  that  vazry  their  size  with  frequency.  The  largest  potential  effect 
exists  at  the  low  end  of  the  frequency  spectrxim  where  the  effective  area  of 
the  antenna  is  very  large  and  can  in  fact  dominate  the  net  Q  of  the  cavity. 

It  is  therefore  better  to  use  probes  or  antennae  which  have  relatively  small 
effective  areas  such  as  B-dot  or  D-dOt  probes. 
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